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We have developed a scheme for a high resolution measurement of the switching current distri- 
bution of a current biased Josephson tunnel junction using a timing technique. The measurement 
setup is implemented such that the digital control and read-out electronics are optically decoupled 
from the analog bias electronics attached to the sample. We have successfully used this technique 
to measure the thermal activation and the macroscopic quantum tunneling of the phase in a small 
Josephson tunnel junction with a high experimental resolution. This technique may be employed 
to characterize current-biased Josephson tunnel junctions for applications in quantum information 
processing. 



I. INTRODUCTION 

The current-biased Josephson tunnel junctioH.|-js an 
ideal system to study both thermal activationnB and 
quantum tunnelingtrH in a controllable experimental en- 
vironment. The dynamics of such a junction is equivalent 
to that of a particle in a tilted washboard potentialtffl. 
The process of the particle escape from a metastable state 
in this system can be characterized, by analyzing the 
transition of the particle from a state in which it is lo- 
calized in a potential well to a state in which it runs 
down the potential. This corresponds to the transition 
of the Josephson junction from the superconducting zero- 
voltage state to a finite voltage state in the presence of 
an applied bias current. At high temperatures the escape 
is dominated by thermal activation across the barrier, at 
low temperatures it is determined by quantum tunneling 
through the barrier. The rate with which the particle es- 
capes from the well depends on the detailed shape of the 
potential, the dissipation in the system and the tempera- 
ture of the thermal bath to which the system is coupled. 

Thermal activation in a current biased Josephson 
junction has been studied botk-. theoretically and ex- 
perimentally—for large dampinguS, intermediate-to low 
dampingEHlErEJ and extremely low dampingE!RL3. At 
low temperatures the quantum mechanical properties of 
Josephson junctions have been investigatad-^Bqrelation 
to macroscopic quantum tuimeling (MQT)ll3'tSt3, energy 
level quantization[XELQ)E3Ej and macroscopic quantum 
coherence (MQC)O. Recently, the quantum mechanical 
properties of Josephson junction systems have regained 
interest in the view of their possible application-iei' solid 
state based quantum information processing' 1 n , l 19 lr c l. The 
prospects of superconducting devices containing Joseph- 
son junctions for use as carriers of quantum information 
has been strengthened by recent encouraging experimen- 
tal resultgB'EaEl. 

In this paper we present a measurement technique im- 
plemented for the characterization of the quantum prop- 
erties of current biased phase qubits. The measurement 
scheme is based on a high resolution measurement of the 
junctions switching currentu using a timing technique. 



The measurement setup was tested by performing mea- 
surements of both thermal activation and quantum tun- 
neling in a current biased small Josephson junction. 

In Sec. II, the Stewart-McCumber model describing 
the dynamics of the phase in a Josephson junction is 
briefly reviewed in relation to thermal activation and 
quantum tunneling of the phase. The measurement tech- 
nique and the setup implemented for its realization is de- 
scribed in Sec. 
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The switching current distributions 
obtained over a large range of temperatures are presented 
in Sec. IV. The relative resolution of these measurements 
is compared with other results from literature in Sec. 
In Sec. VI, the data in the thermal regime is analyzed 



in the view of the low dissipation in the junction. The 
crossover temperature to the quantum regime is deter- 
mined and the measured switching current distributions 
are compared to the quantum pre dictions. Finally the 
results are summarized in Sec. VII. 



II. MODEL 

In the Stewart-McCumber modelEi, the dynamics of 
a current biased Josephson tunnel junction is described 
by an equation of motion for the phase difference <f> 



-I c — sm(4>)-I— - = 0, (1) 
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where $o is the magnetic flux quantum, C is the capaci- 
tance, R the effective junction resistance and I c the fluc- 
tuation free critical current of the junction, see Fig. ^a. 
Here, / is the externally applied bias current. The equa- 
tion of motion (|l|) is equivalent to the damped motion 
of a particle of mass m<j, = C(<I>o/27r) 2 in an external 
potential t/^(</>) along the generalized coordinate (f> 
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In Eq. (Q) the damping coefficient is 1/RC and the po- 
tential is given by 
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FIG. 1: a) Resitively, capacitively shunted junction (RCSJ) 
model, b) U' p versus (j> calculated for a junction with the pa- 
rameters as used in experiments and biased at 7 = 0.995. The 
barrier height Uq and the oscillation frequency lJq are indi- 
cated. Numerically calculated energy levels and the squared 
wave functions are shown. 



where Ej = & n I c /2ir is the Josephson coupling energy 
and 7 = I/I c is the normalized bias current. U^((f)) is a 
cosinusoidal potential (^) with an amplitude proportional 
to Ej, which is tilted proportionally to the applied bias 
current 7. Because of these properties U^((f>) is called 
often washboard potential. 

In the absence of thermal and quantum fluctuation 
and for bias currents 7 < 1, the junction is in the 
zero voltage state, corresponding to the particle being 
localized in one of the potential wells (see Fig. |l|b). 
At finite temperatures T > 0, the particle may escape 
from the well atrJj iag <j aifrr^ nts 7 < 1 by thermally acti- 

by quantum tunnel- 
The rate at which 
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vated processes 
ing through the barrier! 
both processes occur depends on the barrier height 



Uq = 2Ej \fl — 7 2 — 7arccos(7) 
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the oscillation frequency of the particle at the bottom of 
the well 



4 = ^[/"<Vm, = Ml- 7 2 ) 1/4 , 



(5) 



and the damping in the junction. Here uj p = -\/27r/ c /$oC 
is the plasma frequency. For 7 — > 1, Eq. @ can be 
approximated as 



Tjt> 
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At high temperatures the escape of the particle from 
the well is dominated by thermally activated protases, 
which occur with a bias-current dependent rate on 



r t = a t — exp 
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Here a t is a temperature and damping dependent thermal 
prefactor. 



As the temperature T is lowered, the thermal activa- 
tion is exponentially suppressed. At temperatures be- 
low the so called cross-over temperature T* the quan- 
tum tunneling rate of the particle through the barrier 
exceeds the thermal activation rateea. The tunneling 
rate of the phase is calculated in the Wentzel-Kramers- 
Brillouin (WKB) approximation for a particle of mass 
rr$ tunneling through the potential barrier described by 
U^. The effect of damping with the characteristic coeffi- 
cient a = l/(2RCuf) = 1/(2Q), where Q = uj^RC is the 
junction quality factor, is considered it) taons of a cou- 
pling to a bath of harmonic oscillatorsoOCJ. The rate is 
then given by 



T q = Aexp(-S) 
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As in Eq. ([7]), the values of T q , Uq, and ojq depend on /. 

At low temperature and small damping the energy of 
the small oscillations of the phase at the bottom of the 
well is quantized, see Fig. Qb. The energy level quart, 
tization has been , observed experimentally both belowEHI 
and above The escape of the phase from the well 

in the presence of quantized energy levels was discussed 
theoretically by Larkin and Ovchinnikov in Ref. ^2|. This 
theory allows for the calculation of the bias current de- 
pendent escape rate T(I) considering the escape of the 
particle from the well by tunneling from any energy level. 
The occupation of the levels can be calculated for finite 
temperatures and in the presence of microwaves using 
a master equation approach. Previous measurements 
of energy level quantization have been compared with 
the p redic tions of this theory and good agreement was 
found!3E-i This theory has also been used to explain 
switching current distributions for non-stationary distri- 
butions of the phase at high-current ramp rates and above 
the cross-over temperature£d. Here, we analyze our data 
for the escape of the phase in the presence of energy lev- 
els with a separation comparable to the bath temperature 
k^T at temperatures above T* in the framework of this 
theory. 



III. MEASUREMENT TECHNIQUE AND 
SETUP 

The escape of the phase in a Josephson tunnel junction 
is experimentally investigated by performing a statisti- 
cal measurement of the current at which the junction 
switc-hes from the zero-voltage state to a finite voltage 
stateH. There are two different well established meth- 
ods to perform such a measurement. One method con- 
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sists of applying a fixed bias current I < I c to the junc- 
tion and measuring the time between the application of 
the curreat-and the appearance of a voltage across the 
junctionEZru. This corresponds to a direct measurement 
of the lifetime t{I) of the particle in the well. The in- 
verse of the lifetime corresponds to the escape rate r(J). 
In the other method, the bias current applied to the junc- 
tion is ramped up at a constant rate / and the current / 
at which the junction switches .from the zero-voltage to 
a finite voltage state is recordeda. The switching current 
probability distribution P{I) is found by accumulating 
a large number of switching currents I and generating a 
histogram. From the P(I) distribution the bias current 
dependent escape rate 



r(J) = Jin 



(11) 



can be reconstructed^ and compared with the theoretical 
predictions. 



A. Experimental setup 

In order to measure the switching current distribution 
of a Josephson tunnel junction with a high accuracy, we 
have developed a measurement setup in which the ana- 
log biasing electronics and the sample are electrically iso- 
lated from the digital control electronics using an optical 
fiber link. Using a universal time interval counter with 
20 GHz time-base, the switching current is determined by 
measuring the time between the zero crossing of the bias 
current ramped up at a constant rate and the switching 
of the junction to the finite voltage state. 

In the following the elements of the experimental setup 
are discussed starting from the sample cell mounted at 
the cold finger of the dilution refrigerator, the cold and 
warm filters in the dc-bias lines, the room temperature 
analog electronics, the optical fiber link and finally the 
digital control and data acquisition electronics, see Fig. 0. 



1. Sample Cell and Wiring 

The sample, typically fabricated on a 5 x 5 mm 2 chip, is 
mounted in a closed rf-tight copper cell to shield it from 
electromagnetic radiation from the environment. The 
sample cell is thermally anchored to the cold finger of a 
dilution refrigerator. For current-biasing the Josephson 
junction and measuring the voltage across it, four wires 
are fed through microwave filters into the cell. Inside the 
cell the wires are thermally anchored and connected via 
wire bonds to the sample. Each-microwave filter con- 
sists of a commercially availableo 1 meter piece of a 
0.5 mm in diameter, lossy stainless steel coaxial line with 
a 50 51/m inner conductor made from Ni/Cr (80%/20%) 
which is soldered into the wall of [the cell. This type of 
filter, also called thermocoax filtered, has an attenuation 



of more than 50dB/m at frequencies above 1 GHz. Al- 
ternatively, we have also used copper powder microwave 
filtcrst3, consisting of a few meters of resistive wire coiled 
up in a 10 cm long copper tube filled with 10 /im grain- 
size copper powder and mounted in the wall of the cell. 
If heating of the sample at large bias currents is an is- 
sue, the total dc resistance of this filter may easily be 
adjusted by choosing a different wire material and diam- 
eter. This type of filter is designed to have an attenuation 
of more than 50 dB per piece above 1 GHz. All required 
additional wires for heaters, temperature sensors, coils 
etc. were fed into the cell through microwave filters of 
this type. Another possible type of cryogenic microwave 
filter, which has not been tested in our setup, is based on 
a micro-fabricated distributed RLC filtero. 

The sample leads are passed to the 1 K stage of the 
cryostat in a tightly twisted loom of superconducting wire 
and filtered using a low pass RC filtering stage with a 
3 dB cutoff frequency of approximately 50 kHz, which is 
thermally anchored at the 1 K pot. The wiring is further 
fed to the top of the cryostat in a loom of twisted cop- 
per pairs shielded separately from the sensor and control 
wiring of the cryostat. At room temperature all wires are 
going through low pass 7r-type feedthrough filters with a 
cutoff frequency of about 10 MHz. After the last filtering 
stage the wires are separated into shielded pairs which 
are connected to the analog biasing electronics. 



2. Analog Electronics 

The sample is current biased by a voltage controlled 
current source with selectable output-current ranges. 
The source current I is proportional to the input volt- 
age Vj in the range from — 10 V to 10 V. The control 
voltage Vj is a sawtooth waveform with an adjustable 
voltage ramp rate generated by charging a capacitor with 
a constant current. The sawtooth signal is started upon 
a digital trigger signal at the set input (S). The voltage 
ramp is stopped by a digital trigger signal at the reset 
input (R) , see Fig. g. 

The current flowing through the sample is monitored 
by measuring the voltage drop V^Mon across the bias- 
ing resistor Ri using an instrumentation amplifier. The 
voltage V across the junction is measured using a fast 
FET instrumentation amplifier (amplification xlOOO). 
Both voltages Vr,Mcm an d V are fed to two independent 
Schmidt triggers with adjustable threshold and window 
voltages. The current trigger is calibrated to generate a 
TTL pulse Vr ttl upon the zero-crossing of the bias cur- 
rent. The voltage trigger is set up to generate a TTL 
pulse Vttl when the junction switches from the zero 
voltage to a finite voltage state. The trigger window is 
adjusted wide enough to avoid voltage noise induced trig- 
gering. The voltage trigger signal Vttl is simultaneously 
used to stop the current ramp. To avoid a ground loop in 
the signal lines an opto-coupler is used at the reset input 
of the sawtooth generator. 
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FIG. 2: Schematic of the measurement setup. (A detailed explanation is found in the text. 



All analog electronics components described above 
have been specifically designed and implemented for 
these measurements. The supply voltage is delivered by 
lead accumulators. Care has been taken to avoid ground 
loops. All electronics is grounded at a single ground- 
ing point on the top of the cryostat. The cryostat and 
the analog electronics are placed inside a stainless steel 
shielded room. Outside of the shielded room the pump- 
ing lines of the cryostat arc electrically isolated from the 
gas handling system and the pumps. The cryostat is con- 
nected to the shielded room through the pumping lines. 

To perform the current measurement, the current trig- 
ger signal V^ttl and the voltage trigger signal Vttl arc 
converted to digital optical pulses which are passed out of 
the shielded room via optical fibers. Outside the shielded 
room both signals are converted back to TTL level elec- 
tric signals. Conversely, for starting the current ramp a 
TTL pulse is passed via the optical fiber link into the 
shielded room. The arrangement of the analog electron- 
ics setup is shown in Fig. 0. 



3. Digital Control and Data Acquisition Electronics 

A square wave generator provides a TTL signal at a 
rate vi which is used to start the current ramp with the 
same repetition rate z/j. In each cycle of the measure- 
ment the bias current is increased at a fixed rate I until 
the switching current I is reached. The optically decou- 
pled current and voltage trigger signals are then supplied 
to the start and stop inputs of a universal time interval 
counter to measure the time r between the triggers with 
a resolution of approximately 25 ps. Having calibrated 
the current ramp rate / the switching current can be cal- 
culated from the time as I = It. 



Our method has a high instrumental resolution AI = 
J25ps of the current measurement. For typical ramp 
rates between 0.1 A/s and lA/s as used in the experi- 
ments presented here, a relative instrumental current res- 
olution between AI/I = 2.5 1CT 12 and 25. lCT 12 can be 
achieved. This value is much better than that achievable 
using a 16 bit A/D converter for direct current measure- 
ments with AI/I = 0.15 10" 6 . 

Similar time-based methods for measuring the jufttj. 
tions critical current at a constant current ramp rateEflEEl 
or the lifetime of the zero- voltage state of a junction at a 
constant applied bias current! 1 ^rr 5 ! have been implemented 
by other groups. _ 

The time interval counteiO used in our setup allows 
for "the on-the-fly" generation of switching current his- 
tograms from the acquired data during the measurement. 
This feature is very useful for online monitoring of switch- 
ing current histograms with changing experimental con- 
ditions, e.g. when doing spectroscopic measurements 
by applying microwave radiation to the sample or when 
searching for parasitic noise sources. 

The temperature controller of the dilution refrigerator 
is installed inside the shielded room. It can be pow- 
ered either by two lead accumulators or by an isolation 
transformer. The electrical ground of the controller is 
anchored on the cryostat. The temperature controller is 
computer controlled via an optical fiber link. 



IV. EXPERIMENTAL OBSERVATIONS 

The experiments were performed using a high quality 
5x5 /am 2 tunnel junction fabricated on an oxidized silicon 
waver using a standard Nb/Al-AlO^/Nb trilayer process. 
The junction had a nominal critical current density of 
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FIG. 3: Switching current distributions P(I — (/)) at bath 
temperatures between T = 4.2 K and T = 25 mK. The tem- 
perature is color coded according to the scale shown in the 
plot. In the inset the mean value (J) of the switching current 
distribution is plotted versus the bath temperature T. 



j c 1000 A/cm 2 , a nominal capacitance of C « lpF 
and a subgap resistance of R > 500 S7 below 1.0 K. 

The switching current distribution of the sample was 
measured at bath temperatures T between 4.2 K and 
25 mK. The current was ramped up at a constant rate of 
/ = 0.245 A/s with a repetition rate of vj = 500 Hz. Typ- 
ically 10 4 switching currents were recorded at each tem- 
perature. Histograms of the switching currents were cal- 
culated with bin widths of approximately 10 nA in order 
to determine the P(I) distributions. In Fig. || the mea- 
sured switching current distributions are plotted versus 
T. To show all histograms in the same plot, the switching 
probability P has been plotted versus I — (I), where (I) 
is the temperature dependent mean value of the switch- 
ing current shown in the inset of Fig. ||. It is observed 
that the width of the P(I) distributions decreases with 
temperature and then saturates at low temperatures, as 
expected. 

The standard deviation 07 of P(I) is plotted versus 
T in Fig. ^. In the temperature range between 1 K 
and approximately 300 mK, 07 decreases with T, indicat- 
ing the temperature dependent thermal activation of the 
phase across the barrier. At the characteristic tempera- 
ture T* rs 300 mK, 07 saturates at approximately 115 nA 
suggesting that the escape of the phase is dominated by 
quantum tunneling through the barrier. In Sec. VI, the 



P(I) distributions in both the thermal and the quantum 
regime are analyzed quantitatively and compared with 
existing models. 



0.2 0.4 0.6 0.8 1 

bath temperature, T [K] 

FIG. 4: Standard deviation 07 of the P(I) distributions versus 
the bath temperature T. The cross-over temperature T* is 
indicated. 



V. EXPERIMENTAL RESOLUTION 

We have compared the experimental resolution of our 
switching current measurements performed with the ex- 
perimental setup described above to the resolution of 
such measurements presented in literature. Since the 
limiting resolution of the measurement technique is not 
always stated in existing publications, we have focused 
our comparison on the minimum measured width of the 
switching current distribution mentioned in the articles. 
To allow for comparison among very widely varying sam- 
ple parameters, we express the width with respect to the 
critical current of the measured junction. As a measure 
of the width we have considered the standard deviation 
(Tj of the P(I) distributions. If the values of 07 were not 
available in the reference, we evaluated the full width at 
half maximum of P(I). In some of the evaluated ref- 
erences the effect of energy level quantization is investi- 
gated. This typically leads to features on the distribp. 
tions which are more narrow than 07. In such casesE3E£l 
we have taken the width of these smaller features. If 
data of multiple samples were discussed in a reference, 
the data of the sample with the smallest relative P(I) 
width is considered. 

In Tab. | the collected data are summarized. First, we 
note that our sample has a high critical current and a 
relatively small capacitance, leading to a comparatively 
high cross-over temperature. Only in Refs. [i"5|]2"^ higher 
cross-over temperatures have been reported. Most of 
the experiments listed in Tab. | reach measured rela- 
tive widths ai/I c down to a few 10~ 3 . Only in a few 
experiments the smallest relative width is a few times 
10~ 4 . The measurements we present here show a reso- 
lution which is close to the best measurements of this 
type performed so far, see last row of Tab. Q. We stress 
that we have not performed any digital filtering or data 
processing neither on the raw data nor on the switching 
current distributions presented here. In microwave spec- 
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TABLE I: Experimental resolution in Josephson junction switching current experiments. Listed are the reference from which 
data has been extracted (col 1), the focus subject of the respective reference (col 2), the critical current I c of the measured 
sample (col 3), the standard deviation ai of the most narrow measured switching current distribution (col 3), <ji/I c (col 4), the 
junction capacitance C (col 5), the effective junction resistance R (col 6), the cross-over temperature T* (col 7), the junction 
dimensions (col 8), and the materials used for junction fabrication (col 9). 



Refs. 


focus 


Ic 


07 




C* 


R 


T* 


dimensions 


material 






[fxA] 


[ M A] 




[PF] 


[O] 


[mK] \pm x nm] 




Washburn et a7.EZl 


MQT 


57.4 


- 


0.9 x 10" 3 


0.67 


- 


500 


1.2 x 0.13 


Nb-Nb 2 5 -Nb 


Devoret, Martinis et a/HQ 


MQT 


57.4 




1.5 x 10~ 3 


12.00 




150 


1.2 x 0.13 


Nb-Nb 2 5 -Nb 


MQT, ELQ, TA 9.489 




0.4 x 10" 3 


6.35 


190 


37 


10 x 10 


Nb-NbO^-Pbln 


Silvestrini et a/Ej 


TA 


162 


0.245 


1.5 x 10 -3 


500 






100 x 50 


Nb-NbOz-Pb 


Vion et all 


TA 


0.040 




0.4 x 10" 3 


0.00015 








Al-AlO^-Al 


Castellano et aj£j 


TA 


160 




2 x 10~ 3 


0.8 






4x4 


Nb-Al-AlO^-Nb 


Silvestrini et aJtJ 


TA 


47.0 




7 x 10~ 3 


5.5 


20 x 10 3 




10 x 10 


Nb-Al-AKVNb 




ELQ 


47.0 




0.7 x 10" 3 


5.5 


20 x 10 3 




10 x 10 


Nb-Al-AlO^-Nb 


Ruggiero et aJ0 


TA 


175 


0.300 


2 x 10" 3 


2.5 






5x5 


Nb-Al-AlO^-Nb 


this work 


MQT, TA 


315.0 


0.115 


0.4 x 10" 3 


1.00 


500 


300 


5x5 


Nb-Al-AlO^-Nb 



troscopy experiments performed on the same sampleES 
we have observed photon absorption peaks which showed 
even smaller relative widths than the switching current 
distributions presented here. In test experiments we have 
measured current histograms of relative width down to 
1.5 x 1CP 4 , which is the limiting experimental current res- 
olution achievable in our setup at the present time. These 
tests have been performed with the junction biased in 
the resistive state but otherwise with experimental con- 
ditions identical to those used for performing switching 
current measurements. 

We note that the instrumental resolution of the switch- 
ing current measurements by a timing technique as im- 
plemented in this setup is extremely high as in com- 
parison to-|ptker approaches using analog to d igita l 
convertergl3il3^ or time to amplitude convertersE-SEa. 
However, we note that typically the experimental reso- 
lution for switching current measurements is not limited 
by the instrumental resolution of the data acquisition 
hardware but rather by the residual current and voltage 
noise in the junction bias electronics. Thus the most im- 
portant aspect in a well designed measurement setup is 
the reduction of the external interference by the use of 
appropriate filtering techniques together with low noise 
analog electronics. 



VI. DATA ANALYSIS 



the known bias current dependences of the approximated 
barrier height ^) and the small amplitude oscillation fre- 
quency (^|) in the following form 



In 
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(12) 

where we-have introduced the effective escape tempera- 
ture T cs jil The left hand side of Eq. (|l|) is then cal- 
culated in a first iteration using the experimental data 
T(I) and a t ujQ (I) estimated approximating the initially 
unknown value of I c by the maximum measured switch- 
ing current I max and using C = 1.61 pF determined 
from independent microwave .resonance activation mea- 
surements on the same samplec3. Because the right hand 
side of Eq. (|l2|) is a linear function of the bias current /, 
both the critical current I c in absence of fluctuations and 
the effective temperature T esc of the thermal activation 
process can be determined from a linear fit of the exper- 
imental data to the right hand side of Eq. (|l2|). Accord- 
ingly, we find 



T 



Cconst 
Clin 

1 $ 4\/2 



1 



A. Thermal Activation 



k B 2tt 3 1/2 n . 



(13) 
(14) 



At temperatures T > T* the escape of the phase is 
dominated by thermal activation. At a fixed bias current 
the activation rate is given by Eq. (^). Using Eq. (|Tl"|), 
the experimental escape rate T(I) is determined from the 
P{I) data and compared with the predictions of Eq. (]?]). 

To perform this comparison we rewrite Eq. (uh using 



where c CO nst and cn n are the two fitting parameters. The 
values of I c and T osc are then found with high accuracy by 
iteratively repeating the fitting procedure with the value 
of I c found in the previous iterationE£l. This procedure is 
converging quickly due to the logarithmic dependence of 
the left hand side of Eq. ( |l2| ) on I c . 
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1. The thermal pref actor in the extremely low damping 
limit 

We have extracted the effective escape temperature 
T csc from the experimental data according to the pro- 
cedure described above taking into account the thermal 
prefactor at- In the simple transition state theory, the 
thermal prefactor is unity and the escape rate is only de- 
termined by the fraction of particles with energy larger 
than the barrier height. In the high as well as in the 
low dampiug regimes the prefactor deviates substantially 
from unitycJ. 

The exact form of the prefactor is found from a analysis 
of a classical Langevin equation for the motion of the par- 
ticle according to Eq. (g) in the presence of fluctuations 
(see e.g. Refs. p|p^). In the moderate to strong damping 
limit the value of at depends predominantly on the damp- 
ing and the small oscillation frequency in the junction. 
In the very strong damping limit, the inertia of the par 
cle can be neglected and its motion becomes diffusiveL 
In the low to extremely low damping regime, however, 
the distribution of the particles in the well deviates from 
thermal equilibrium. The thermal escape across the bar- 
rier depletes the population within a range of k^T to the 
top of the barrier. Therefore, at becomes dependent on 
the barrier height and on temperature. For simplicity, we 
use an approximate analytical expression for the thermal 
prefactor in the extremely low to low damping regime, 
relevant for the exp erimental data presented here, as cal- 
culated in Ref. 40 



a i 



4dn 



aoQk B T 



1.8 W 



(15) 



where ao is a numerical constant close to unity. An ex- 
act calculation of at valid in all damping regimes can be 
found in Ref. § 

In a first step, T esc has been evaluated for at = 1. Un- 
der this assumption the escape temperatures extracted 
from the data deviates more than 100 mK from the bath 
temperature at T > T* , see open symbols in Fig. [|. The 
difference T esc — T even increases with increasing tem- 
perature, as shown in the inset of Fig. ||. This deviation 
is too large to be due to an inaccurate estimate of the 
junction capacitance or errors in the calibration of the 
current measurement, thus, indicating the importance of 
the thermal prefactor. 

In a second step, we used the thermal prefactor given 
by the expression ({f|) in the calculation of the left hand 
side of Eq. ( |l2| ) in our data analysis scheme. We used 
the junction resistance R as the adjustable parameter 
determining the junction quality factor Q. We found 
good agreement between the extracted T osc and the bath 
temperature T for R = 500 f2 neglecting its tempera- 
ture dependence, see solid symbols in Fig. 0. With the 
low damping prefactor (|i"5|), the difference \T CSC — T\ is 
less than 75 mK in the temperature range above T* , as 
shown by the solid symbols in the inset of Fig. pi The 
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FIG. 5: Escape temperature T osc extracted from experimental 
data versus bath temperature T for at — 1 (open symbols) 
and at according to Eq. ( [l5| ) (solid symbols). At high T the 
statistical error in T esc extracted from the fitting procedure is 
approximately given by the symbol size, at low T it is substan- 
tially less than the symbol size. The solid line corresponds to 
T = T csc . The inset shows T esc - T. 
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effective resistance used for the fit is in good agreement 
with the sub-gap resistance determined from the dc cur- 
rent voltage characteristic of the junction. 

In Fig. ||, the thermal prefactor at evaluated at the 
most probable switching current I p is plotted versus tem- 
perature for T > T* . Independently of the bath temper- 
ature, the mean value of a t is approximately 0.17 . This 
weak temperature dependence can be attributed to the 
almost constant ratio of k^T /Uq(I p ) for the switching 
data measured at different temperatures, see Eq. JTe|). 
At a fixed bath temperature, the bias current variation 
within a single P(I) distribution results in a noticeable 
variation of at- The value of a t at the minimum and 
maximum switching currents J m i n and / max is indicated 
in Fig. ^ by error bars. 
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Thus, the experimental data is explained with good 
accuracy by the classical model of the thermal activation 
of the phase over a potential barrier in the presence of 
extremely small damping. The quality factor used for 
the extraction of the effective escape temperature is in 
agreement with the value determined from the sub gap 
resistance of the junction. 



2. Thermal activation in the presence of energy levels 

In thermal activation theory it is assumed that the 
particle energies in the well are continuously distributed. 
For the sample considered here, however, the quantiza- 
tion of the energy of the small amplitude oscillations of 
the phase may be relevant. In fact, we found direct ev- 
idence for the energy level quantization of the phase in 
this sample by using microwave spectroscopy^. 

To address this issue we calculate the approximate 
number of levels in the well at a given bias current as 



20 f 



(16) 



TJ<P 

Tiuji 



neglecting the anharmonicity of the well. Eq. (|16|) has 
been evaluated at the most probable switching current 
I p at each temperature, see Fig. f^a. In the quantum tun- 
neling regime at T < T* , we find that there are only a few 
levels in the well (n ma x < 3). The number of levels in the 
well at the relevant switching currents increases approx- 
imately linearly with temperature from about 5 levels 
at IK to about 20 levels at 4K. Additionally we have 
evaluated the level separation AE — JiuJq(I) in harmonic 
approximation at the most probable switching current I p , 
see Fig. |^b. At temperatures below 1 K the level separa- 
tion is much larger than k^,T and at temperatures above 
1 K it is comparable with k^T. These facts indicate that 
the energy level structure may be relevant for the process 
of the escape of the phase even at temperatures T > T* . 

To analyze the effect of the level quantization on the 
escape of the phase at high temperatures we have made 
use of the Larkin Ovchinnikov theoryH. In the frame- 
work of their model the escape of the phase is considered 
in terms of tunneling out of the well from individual en- 
ergy levels. The overall escape rate Tlo(I) depends on 
the population of different levels, which is determined by 
the coupling to the thermal bath. To apply the Larkin 
Ovchinnikov theory, we have calculated the energy level 
structure of the phase in the potential in dependence on 
the bias current /. The matrix elements for transitions 
between individual levels have been evaluated and the 
resulting escape rate has been calculated by solving a 
master equation for the dynamics of the system. From 
the escape rate Flo (I) we have calculated P(I) distri- 
butions in a wide temperature range. 

The calculated P(I) distributions at temperatures well 
below T* , when essentially only the ground state is popu- 
lated, were in good agreement with the distributions cal- 
culated using the WKB approximation for the tunneling 
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FIG. 7: a) Approximate number of quantized energy levels 
n max in the well at the most probable switching current I p 
versus the bath temperature T. b) Harmonic approximation 
of energy level separation AE/Icb at I p in Kelvin versus T. 



from the ground state. At temperatures above T*, the 
width of the distribution calculated using Larkin Ovchin- 
nikov theory increased with temperature as expected. We 
have evaluated these calculated data in the same way as 

We found that 



the experimental data, see Sec VI A 1 



the calculated histograms are consistent with the classi- 
cal prediction for a low damping thermal prefactor. Thus 
in the high temperature limit the results of the Larkin 
Ovchinnikov theory are in agreement with the predic- 
tions of purely classical thermal activation theory. .More 
details on this analysis will be presented elsewheraid. 



B. Quantum Tunneling 

At temperatures below the crossover temperaturcS 



27rfcr 



1 



> 2 \ 1/2 
1 \ \ 1 



2Q, 



2Q 



(17) 



the escape of the phase from the potential well is dom- 
inated by quantum tunneling through the barrier. Ac- 
cording to Eq. (|l7|), the crossover temperature T* for 
this sample at the most probable switching current I p is 
predicted to be approximately 282 mK, which is in good 
agreement with our experimental findings, see Fig. |J. 
Since the quality factor of this sample is bounded by 
Q > 100, the reduction of T* in comparison to the ideal 
case of Q = oo is less than one percent. Accordingly, 
the reduction of the tunneling rate due to dissipation is 
small. 

Furthermore, we have compared the measured P(I) 
distributions for T <C T* to the predicted distributions 
calculated using the quantum tunneling rate (||). As an 
example, the data for T = 25 mK are shown in Fig. |8| 
by solid symbols. The quantum tunneling distribution 
P q (I) calculated for I c = 325.05 fiA and C = 1.61 pF 
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FIG. 8: Measured switching current distribution P(I) at 
T = 25mK (solid symbols). The solid line is the predicted 
switching current distribution P q (I) in the quantum tunnel- 
ing regime for I c = 325.05 /xA and C = 1.61 pF. The P q (I) 
distribution broadened by residual experimental noise mod- 
eled by gaussian distribution with ai = 47.5 nA (dotted line) 
is shown by the dashed line. 

is shown by the solid line. The critical current I c was 
extracted from fits to the data in the thermal regime. 
As it has been already said above, the capacitance C 
was determined with high accuracy in independent spec- 
troscopic measurjements of the energy level structure of 
the same sampleEl The most probable switching cur- 
rent of the calculated distribution is in good agreement 
with experimental data. However, the width of the mea- 
sured distribution is slightly bigger than the predicted 
one, which is most likely due to the finite residual elec- 
tromagnetic interference in the measurement setup. To 
quantify this effect, we have calculated the convolution of 
the predicted quantum tunneling distribution P q (J) with 
a gaussian current distribution with a standard deviation 
of (7/ = 47.5 nA, shown by the dotted line centered about 
I = 323.4 /zA. The resulting distribution, shown by the 
dashed line in Fig. |§|, accurately fits to the experimental 
data. 



VII. CONCLUSIONS 

We have developed an experimental setup for the in- 
vestigation of the escape of the phase in a current-biased 



Josephson tunnel junction. Using a timing technique 
we are able to achieve a high instrumental resolution in 
the switching current measurement of a Josephson junc- 
tion. The relative resolution of the switching current 
measurements demonstrated in this paper is better than 
3.7 x 1CP 4 . This distribution width aj/I c is comparable 
to the highest resolution measurements of this type per- 
formed so far. Using the described technique we have in- 
vestigated the escape of phase in a high-quality 5x5 /im 2 
Nb-Al/AlO^-Nb tunnel junction with a relatively large 
critical current of I c rs 325 /iA. At temperatures above 
the crossover temperature of T* » 300 mK the escape is 
dominated by thermal activation. Due to the high qual- 
ity factor of the junction (Q > 100) the thermal activa- 
tion is observed in the limit of extremely low damping. 
The thermal prefactor at extracted from the experimen- 
tal data is varying between 0.1 to 0.3 depending on bias 
current and temperature. At temperatures below T*, 
the value of which agrees well with the theoretical pre- 
dictions, the escape of the phase by quantum mechanical 
tunneling is observed. Due to the high Q-factor of the 
junction, the tunneling rate is only weakly suppressed by 
dissipation. 

Recent developments in the field of quantum informa- 
tion processing have strongly increased the interest in 
research on the cpantum properties of Josephson junc- 
tion systemsEJoH Using the measurement technique 
and setup described in this paper the quantum mechan- 
ical properties of current biased Josephson junctions can 
be investigated with high accuracy. In particular, we are 
currently applying this technique for the experimental in- 
vestigation of the theoretically predictedc3c3 but not yet 
observed quantum properties of a single vortex in a long 
Josephson junctions. 
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